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Abstract 


Certain  spread-spectrum  systems  make  use  of  waveforms  built  of  basic 
pulses  spread  by  PN  modulation  of  some  kind.  In  this  study  the  two  most 
common  modulation  waveforms,  BPSK  and  MSK  are  discussed  in  detail.  The 
objective  is  to  derive  receiver  structures  utilizing  sampled-data  inte¬ 
gration  of  the  chips,  for  example,  using  CCD  or  digital  correlators.  A 
constructive  design  procedure  is  described,  together  with  an  analysis  of 
performance  in  terms  of  SNR  loss  due  to  sampling.  It  is  shown  that  sampling 
loss  can  be  reduced  by  proper  choice  of  the  prefilter  which  precedes  the 
sampling  of  the  data.  _ 
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Summary 

This  report  is  concerned  with  certain  aspects  of  the  design  of  receivers 
for  two  waveforms  commonly  used  as  spreading  modulations  in  jam-resistant 
communications  systems.  Specifically,  we  are  interested  in  the  reception  of 
bursts,  each  many  chips  long,  in  which  a  PN  code  controls  the  modulation  by 
either  phase  changes  of  180°  (BPSK)  or  carrier  frequency  shifts  as  in  MSK. 
Complex  signal  structures  can  be  built  on  these  bursts.  In  all  cases  the 
signal  parameters  (including  the  code)  are  known  to  the  receiver,  and  the 
problem  of  interest  is  to  optimize  its  design. 

It  is  natural  to  base  the  receiver  design  on  the  matched  filter 
principle,  insofar  as  this  is  possible  with  the  hardware  available,  and  one 
of  the  features  of  this  study  is  a  description  of  the  form  taken  by  a  matched 
filter  for  waveforms  of  this  type.  The  optimum  receiver,  as  it  might  be 
realized  at  base  band,  forms  in-phase  (I)  and  quadrature  (Q)  baseband  signals 
and  passes  each  through  a  filter  matched  to  a  pulse  shape  characteristic  of 
the  waveform.  This  pulse  is  a  rectangle,  one  chip  long,  for  BPSK,  and 
one-half  cycle  of  a  sinusoid,  two  chips  long  overall,  for  MSK.  These  outputs 
are  then  passed  through  delay  structures  which  match  the  chip  structure  of  the 
waveform.  Each  chip  delay  is  tapped  and  weighted  with  the  corresponding  code 
bit  and  summed.  The  final  I  and  Q  outputs  are  summed  and  thresholded  for 
sync  detection  or  data  demodulation. 

The  optimum  delay  structure  is  continuous  in  nature,  such  as  a  delay 
line.  Many  systems  use  sampled  data,  however,  either  analog  or  quantized  in 
amplitude,  and  the  sampled  base-band  filter  outputs  are  then  passed  through  a 
discrete  delay  structure,  such  as  a  CCD  correlator  or  a  digital  correlator,  as 
a  practical  approximation  to  the  optimum  structure. 

In  sampled-data  systems  of  this  kind,  performance  varies  with  "sampling 
phase",  the  relative  timing  of  the  sampling  device  and  the  chip  boundaries  of 
the  actual  received  waveform.  Pure  matched  filter  performance  is  attained 
only  for  certain  values  of  the  sampling  phase,  and  on  the  average  a  certain 
loss  in  performance  occurs,  usually  expressed  as  a  "sampling  loss”  in 
effective  signal-to-noice  ratio  (SNR). 

The  sampling  loss  is  the  main  concern  of  this  study,  and  it  is  evaluated 
as  a  function  of  sampling  rate  (number  of  samples  per  chip)  for  both 
waveforms,  using  the  baseband  filters  matched  to  the  appropriate  pulse,  as 
described  above.  However,  it  is  also  shown  that  this  sampling  loss  behavior 


vi 


can  be  changed,  often  in  a  desirable  way,  by  modifying  the  baseband  filters, 
and  several  filter  designs  are  presented,  based  on  different  criteria  for 
acceptable  sampling  loss  behavior. 

The  reductions  in  sampling  loss  are  not  great,  but  it  is  shown  that  the 
sampling  loss  variation  with  sampling  phase  can  be  sharply  reduced  for  BPSK, 
even  with  one  sample  pe.’  chip,  and  that  sampling  losses  with  HSK  sampled  once 
per  chip  can  be  comparable  to  losses  with  BPSK  sampled  at  twice  the  rate. 

These  conclusions  should  be  useful  in  receiver  design,  even  though  they  deal 
with  only  one  aspect  of  that  design,  since  sampling  rates  (and  chipping  rates) 
are  usually  pushed  as  high  as  the  technology  will  allow. 

In  addition  to  these  details,  a  number  of  mathematical  properties  of  MSK 
are  discussed,  including  the  desirability  of  using  one  of  the  two  tones  as  a 
reference  for  demodulation  to  baseband. 


I.  introduction 

Many  communications  systems  make  use  of  waveforms  built  entirely  of  basic 
pulses,  combined  in  ways  to  achieve  the  fundamental  requirements  of  synchro¬ 
nization  and  information  transmission.  Synchronization  implies  accurate 
timing  of  the  arrival  of  a  reference  point  in  the  waveform,  in  support  of 
data  demodulation  and,  in  some  cases,  ranging.  If  data  modulation  takes  the 
form  of  pulse  position  modulation,  then  the  entire  task  of  the  receiver  reduces, 
in  a  sense,  to  the  detection  and  timing  of  the  arrival  of  the  individual  basic 
pulses  of  the  waveform. 

One  approach  to  waveform  design  is  to  use  intense,  very  short  pulses,  the 
other  common  method  is  to  increase  the  bandwidth  of  a  pulse  of  unmodulated 
carrier  by  modulating  it.  Some  form  of  phase  modulation  is  often  preferred 
because  of  the  convenience  of  dealing  with  constant-envelope  pulses. 

Binary  phase  modulation  schemes  have  a  natural  advantage  in  that  the  wave¬ 
form  itself,  as  well  as  the  matched  filter  or  correlator  included  in  the  pre¬ 

ferred  implementation  of  the  receiver,  can  be  neatly  specified  by  means  of  a 
"code",  or  sequence  of  bits.  In  the  applications  of  this  approach  to  waveform 
design,  the  binary  code  is  used  only  to  spread  the  pulse  spectrum  and  not  to 
convey  information,  since  the  code  is  already  known  to  the  receiver.  Informa¬ 
tion  could  be  transferred  by  transmitting  one  of  a  set  of  codes,  but  this 
forces  the  receiver  to  implement  parallel  channels,  each  designed  to  detect 
one  of  the  codes,  and  hence  each  channel  is  basically  similar  to  a  fixed-code 
receiver,  with  analogous  requirements  on  the  pulses. 

Two  of  the  simplest  modulation  schemes  which  translate  these  ideas  and 
constraints  to  practice  are  binary  phase  shift  keying  (BPSK)  and  minimum- 
shift  keying  (MSK) .  In  each  case  the  pulse  is  divided  into  a  number  of  equal 
segments,  called  chips,  one  for  each  code  bit.  In  BPSK,  the  carrier  phase 

takes  one  of  two  values  in  each  chip,  according  to  the  corresponding  code  bit, 

while  in  MSK  the  instantaneous  frequency  takes  one  of  two  values,  in  analo¬ 
gous  fashion.  In  BPSK,  the  two  values  of  phase  differ  by  180°,  while  in 
MSK  the  difference  of  the  two  frequencies  accumulates  a  phase  difference  of 
180°  over  the  duration  of  one  chip. 
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This  report  is  a  study  of  a  class  of  receiver  structures  designed  to  de¬ 
tect  and  time  basic  pulses,  spread  by  either  BPSK  or  MSK.  The  chief  charac¬ 
teristic  of  these  receivers  is  their  use  of  sampled-data  correlators  for  the 
implementation  of  a  portion  of  the  processing,  and  a  main  objective  of  the 
study  is  the  optimization  of  the  remaining  portion  of  the  receiver  in  order  to 
minimize  the  inherent  loss  in  performance  due  to  sampling.  The  emphasis  is  on 
detection  performance,  but  the  results  obtained  are  also  useful  in  assessing 
timing  accuracy  as  well.  In  addition,  some  insight  is  gained  into  system 
trade-off  issues,  relating  sampling  rate  and  correlator  structure  to  perfor¬ 
mance. 


II.  THE  MATCHED  FILTER 

In  this  section  we  review  the  simple  theory  of  the  matched  filter. 

This  will  serve  to  establish  some  basic  results  in  the  general  case,  avoid¬ 
ing  repetition  for  specific  waveforms,  and  it  will  also  introduce  some  con¬ 
ventions  of  notation. 

The  received  waveform  is  represented  as  a  modulated  carrier: 

R(t)cos[coot  +  $(t)]  , 

and  the  function  we  deal  with  is  the  complex  modulation 
Z(t)  =  R(t)e1$(t) 

This  representation  is  not  unique,  of  course,  since  the  same  waveform  can  be 
expressed  in  terms  of  a  different  carrier  frequency,  as  follows: 

R(t)cos  [co^t  +  $'(t)]  , 

where 

$'(t)  =  *(t)  +  (w  -  w'  )t  • 
o  o 

This  is  a  simple  point,  of  course,  but  in  the  discussion  of  MSK  it  proves 
useful  to  consider  two  such  representations,  with  slightly  different  carrier 
frequencies . 

The  real  and  imaginary  components  of  Z(t)  are  recovered,  as  in-phase  (I) 
and  quadrature  (Q)  signals,  by  beating  the  received  signal  to  baseband  with 
a  local  oscillator  at  the  chosen  carrier  frequency  (together  with  a  90-degree 
phase  shifter  to  provide  the  quadrature  reference). 

The  complex  modulation  is  represented  as  a  sum  of  signal  plus  noise: 

Z(t)  =  Ael8  S(t-T)  +  N(t)  . 

Here,  A  and  B  represent  amplitude  and  carrier  phase,  while  S(t)  describes 
the  basic  signal  waveform,  usually  normalized  to  unit  amplitude.  The  para¬ 
meter  T  represents  the  actual  signal  arrival  time,  and  N(t)  describes  the 
complex  noise.  This  noise  is  assumed  to  be  white  with  the  properties 
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E  N(t)  =  E  N(t)  N(t')  =  0  , 

E  N*(t)  N(t')  =  2Nq  6(t'-t) . 

The  real  and  imaginary  parts  of  N(t)  are  uncorrelated,  and  N  is  the 

o 

single-sided  power  spectral  density.  E  stands  for  expectation  value,  or 
ensemble  average. 

The  output  of  the  matched  filter  to  an  input  Z(t)  is  taken  to  be 
W(t)  =  Js*(  o  )Z(  o  +  t)da  . 

For  convenience,  we  have  expressed  the  matched  filter  output  as  a  complex 
correlation  with  the  signal  waveform,  rather  than  a  convolution  with  the 
corresponding  filter  impulse  response.  Also,  no  attempt  is  made  to  account 
for  realizability  delays.  We  also  use  the  convention  that  all  integrals 
run  from  -°°  to  +  °°,  since  in  every  case  the  integrals  contain  signals  of 
finite  duration  as  factors. 

We  are  concerned  here  with  incoherent  reception,  hence  detection  and 
timing  of  signals  are  based  upon  the  magnitude,  |w(t)j,  obtained  as  the  root- 
sum-square  of  the  I  and  Q  components  of  the  filter  output.  If  the  same  re¬ 
ceived  waveform  is  represented  in  terms  of  another  carrier  frequency,  to'o, 
the  modulation  will  be 

Z'(t)  =  el(uo_a)'o)t:  Z(t)  , 

and  the  signal  component  will  be 

S'(t)  =  ei(wo-UJ'o)t:  S(t)  . 

The  corresponding  matched  filter  output  will  be 

W’(t)  =  /s'*(o  )  Z'(o  +  t)do  =  el(u,o-a)’o)t:  W(t)  . 

Since  W' (t)  has  the  same  magnitude  as  W(t),  there  will  be  no  difference  in 
performance  of  the  filters  corresponding  to  the  two  choices  of  carrier,  and 
again  this  will  be  useful  in  analysis  of  receivers  for  MSK. 


4 


In  the  absence  of  noise,  the  onto”*-  „ouid  be  due  to  signal  only: 


W(t) 


Ae1^/s*( o  )S( a  +  t  -  x)d o 
Ae^  C(t  -  t)  , 


where  C(t)  is  the  signal  autocorrelation  function: 

C(t)  =  fs  (o)S(o+  t)do  . 

This  output  has  its  peak  value  at  t=T: 

W(t)  =  Ael6C(0)  . 

For  noise  alone,  the  output  is 

W(t)  =  /s*(c)N(o  +  t)d  o  , 
a  random  process  with  autocorrelation  function 


EW*(t)  W(t')  =  2Nq  jJs(o  )S*(  d)8(a  +  t-o'-t')do  d  o  ' 

=  2N  C(t'-t)  • 
o 

i  1 2 

The  output  signal-to-noise  ratio  (SNR)  is  defined  as  the  ratio  of  |W(t)|  , 

2 

for  signal  only,  to  E | W( t ) [  ,  for  noise  alone.  The  peak  SNR,  which  occurs 
when  t=T,  is 

St™  -  -  VV 

o 

where  Eg,  the  total  signal  energy,  is 


Es  =  J-  c(0)  =  J-  JjS(o  )  |  2  do. 


The  signals  of  interest  in  this  report  are  of  constant  envelope:  [ S ( t )  |  =  1, 

2 

and  of  finite  duration,  say  T,  so  that  Eg  =  (1/2)  A  T. 

In  later  sections  we  will  be  discussing  approximate  implementations  of 
matched  filters,  and  their  performance  will  be  characterized  by  the  loss  in 
peak  output  SNR,  relative  to  E^/N^  . 
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III.  THE  BPSK  SIGNAL  AND  MATCHED  FILTER 


The  BPSK  signal  can  be  described  by  the  modulation  function 


S(t)  =  E  a  P  (t-nA) 
n  o 
n 


The  a^  are  binary  variables,  assuming  values  ±  1,  for  n  in  the  range  1  <  n  <  L. 

They  represent  some  pseudorandom  code,  whose  detailed  properties  are  not  of 

interest  here,  except  for  sidelobes,  as  discussed  below.  We  use  the  convention 

that  a  =0  for  n  <  0  and  n  >  L,  so  that  all  sums  are  unrestricted  in  extent, 
n 

The  chip  duration  is  A,  and  the  basic  pulse,  Pq  is  defined  by 


P  (t) 
o 


—  <  t  <  ^ 
2  ~  '  2 


otherwise 


It  is  unnecessary  to  make  any  assumption  about  the  relation  between  carrier 
frequency  and  chip  duration,  since  the  analysis  uses  only  the  modulation. 

In  practical  schemes  for  generation  of  BPSK  it  may  be  useful  to  work  at  a 
carrier  which  is  a  multiple  of  the  chipping  rate,  but  this  relation  is  lost 
when  the  carrier  is  changed,  either  before  transmission  or  in  the  receiver 
as  part  of  the  demodulation  process. 

In  terms  of  p  ,  the  normalized  autocorrelation  function  of  P  : 

0  o 

po(t)  =  1  IVo(o  }  P0(0+  t)d0 

/  I  I 

I1  -  |  ;  -A  <  t  <  A 

0  ;  otherwise  , 


the  BPSK  signal  autocorrelation  function  is 
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C ( t )  =  E  anam  fPo(  0  -  n  A) Pq ( a  +  t  -  m  A)do 

n,m  J  ° 


=  A  E  C£  pQ(t  -  £A)  . 

We  have  introduced  the  code  autocorrelation  sequence 


=  E  a  a 

n  n+£ 
n 


Note  that  Cq  =  L,  C  ^  and  that,  for  Z  >  0, 

L-£ 

For  "good"  code  sequences,  the  sidelobe  values 
(C^  for  £  ^  0)  will  be  small  compared  to  L.  In  any  case. 


C.  =  E  a  a  ,  „ 
£  ,  n  n+£ 

n=l 


so  that  =  0  for  |£|_>  L. 


C(0)  =  A  Cq  pQ(0)  =  L  A  , 


which  equals  to  signal  duration,  hence 


E  =  7t  A2  L  A  . 
s  2 


The  matched  filter  processor  takes  a  simple  form,  since  we  may  write 


where 


W(t)  =  E  a^  J  PQ (  o  -  n  A)  Z  ( a  +  t)do 
n  J 

=  E  a  Z  (t  +  n  A)  , 
n  i 
n 


Z,  (t)  =  f P  (a)  Z(a  +  t)da  . 
1  ./  o 


7 


Now  (t)  can  be  obtained  by  passing  Z(t)  through  a  filter,  in  fact  a  filter 
matched  to  the  pulse  P  ,  again  with  the  realizability  delay  ignored.  It  is 
immaterial  whether  this  filtering  is  accomplished  at  bandpass,  with  the  signal 
on  a  carrier,  or  as  a  pair  of  identical  filters  (since  Pq  is  real)  operating 
on  the  baseband  I  and  Q  components  of  Z(t).  The  desired  output,  W(t),  is  ob¬ 
tained  from  the  prefiltered  waveform,  Z^(t),  by  passing  the  latter  through 
a  tapped  delay  line  structure.  Each  segment  of  this  line  introduces  delay  A, 
and  the  tapped  outputs  are  weighted  by  the  code  weights,  a^,  and  then  summed. 
Again,  our  formula  expresses  the  output  as  a  correlation  instead  of  a  convolu¬ 
tions,  and  the  realizability  delay  is  ignored.  The  delay  line  structure  can 
be  implemented  on  a  carrier  or  (since  the  a^  are  real)  as  a  pair  of  identical 
structures  operating  on  the  prefiltered  baseband  components  of  Z^(t). 

The  output  SNR  of  the  matched  filter,  as  a  function  time,  will  be 


SNR(t) 


A2  | C (t-T) | 2 

2N  C(0) 
o 


%  C(t-T) 
Nq  C(0) 


E  C.  p  (t-x-£A) 
i  *  0 


The  expression  inside  the  bracket  takes  the  value  unity  at  t  =  x,  the  signal 
arrival  time,  and  the  values  C^/L  at  times  t  =  T  +  £A.  Moreover,  this  sum 
varies  linearly  between  these  values  at  intermediate  times,  as  is  easily  seen 
by  substituting  for  pQ(t).  Since  the  filter  produces  a  continuous  output  in 
time,  the  peak  value,  Eg/No,  of  SNR(t)  is  always  attained. 
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IV.  THE  MSK  SIGNAL  AND  MATCHED  FILTER 


We  introduce  the  MSK  modulation  as  pure  frequency  modulation,  where 
within  each  chip  the  radian  frequency  is  either  increased  or  decreased  by  a 
fixed  amount,  v,  according  to  the  value  of  a  code  bit.  The  value  of  v  is 
related  to  the  chipping  rate  so  that  the  modulation  produces  a  phase  change 
of  ±  90  degrees  over  one  chip.  We  put 

1i$( t)  ;  0  <  t  <  L  A 

‘ 

0  ;  otherwise 

where  the  initial  phase,  $(0) ,  is  zero,  and  the  instantaneous  frequency  is 


<f(t)  =  bn  V  ’  ^or  (n--^  A  £  t  <  n  A  . 

The  are  binary  variables  representing  a  code  sequence,  and  n  runs  from 
1  through  L.  The  relation  between  modulation  frequency  and  chip  length  is 
simply 


v  A  =  tt/2  . 

The  phase,  4>(t),  is  a  continuous  function  of  time,  which  takes  on  the 
values  $(A)  =  b^TT/2  ,  $(2A)  =  (b ^+b /2 ,  and  generally 

$(nA)  =  (b.  +  ...  +  b  )tt/2  , 

i  n 

at  the  chip  transition  times.  At  an  intermediate  time  the  phase  varies 
linearly: 


$(t)  =  $(nA)  +  b  ,  v(t-nA)  , 

IVri 


for 


n  A  <  t  <  (n+l)A  . 


For  a  binary  variable,  b,  we  have  the  obvious  identity 


b  ®  =  cos0  +  i  b  sin0  , 


and  in  particular 

i  b  */2  <  v 

e  “  1  b 

Therefore,  the  modulation  function  can  be  written  in  the  form 

S(t)  =  e1  (bl+  •••  +  bn)(7T/2)  +  1  bn+lv(t  "  nA) 

=  inb^ • • *bn  [cos  v  (t  -  nA)  +  i  bn+1  sin  v(t  -  nA)], 

for  t  in  the  range  [n  A,  (n+l)A].  This  formula  holds  for  t  in  the  range 
[0,  A]  if  we  interpret  the  product  b^...b  to  be  unity  for  n  =  0.  We  define 
the  new  sequences  of  binary  variables: 


a  =1 
o 


a  =  b,  . . .b  ,  l<n<L  , 

n  In’  —  —  ’ 


and  observe  that 


S(t)  =  inancos  V  (t-nA)  +  in+^  a^+^  cos  V  [t  ~(n+l)A]  , 

for  t  in  the  range  [n  A,  (n+l)A]  and  0  n  <  L.  We  have  used  the  relation 
of  v  to  A  to  equate  sin  V  t  to  cos  V(t-A)  . 


If  we  define  the  "MSK  puise",  P(t),  by 

'  cos  vt  ;  -A  _<  t  _<  A 

;  otherwise  , 


P(t)  = 


!cc 
0 
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t  hen 


P[t-(n  +  1)A]  , 


S(t)  -  ina  P(t-nA)  +  in+1a 

n  n+1 

for  t  In  [nA,(n  +  1 ) A] .  But  *.or  t  in  this  range,  the  sum  can  be  extended: 

L 

S(t)  -  £  ina  P(t-nA)  , 
n 

n»o 

since  the  other  terms  vanish  by  the  definition  of  P(t).  Obviously,  this  sum 
is  a  correct  expression  for  S(t)  for  any  t  in  whole  range  [0,LA],  and  it  serves 
to  represent  the  MSK  modulation  as  a  sum  of  overlapping  pulses  with  complex 
weights  of  a  particular  kind.  It  should  be  noted  that  the  code  sequence  used 
in  this  representation  is  not  the  original  sequence,  which  is  easily  recovered 
from  the  identity 

bn  =  Vn-l  *  (1  <  n  <  L)  . 

This  representation  is  so  useful,  both  for  analysis  and  as  a  model  for  the 
■'ctual  generation  of  MSK,  that  we  redefine  S(t)  to  be  the  value  of  this  sum  for 
all  t,  even  outside  the  interval  [0,LA].  This  means  that  we  allow  the  wave¬ 
form  to  depart  from  the  constant-envelope  form  in  the  end  intervals  [-A,  0] 
and  [LA,(L  +  1)A],  where  a  single  cosine-term  spills  over.  We  take  the 
a^-sequence  to  be  the  basic  code,  and  make  one  final  change,  by  dropping  the 
term  n  *  o,  so  that 

L 

S(t)  =  £  ina  P(t-nA)  , 

t  n 

n*l 

which  vanishes  outside  the  range  [0,(L  +  1 ) A ] .  The  envelope  is  unity  within 
[A,  LA],  and  it  is  easy  to  see  that 

J  [S(t)|2dt  =  LA  . 

By  making  this  last  change  we  gain  a  consistency  with  the  BPSK  formulation, 

and  we  also  use  the  convention  that  a  =0  unless  1  <  n  <  L,  so  that  sums 

n  —  — 

can  be  unrestricted  in  range. 
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If  the  carrier  radian  frequency  is  a)o,  then  the  MSK  waveform  has  an 
instantaneous  frequency  of  either  u>q  +  V  or  -v.  Looked  at  another  way, 
the  actual  waveform  presents  one  of  two  frequencies,  which  differ  by  half 
the  chipping  rate,  and  our  representations  of  the  modulation  refer  to  a 
carrier  which  is  the  average  of  these  two  frequencies.  Using  the  freedom 
to  change  the  carrier  frequency  to  represent  the  same  real  waveform,  which 
was  discussed  in  Section  II,  we  can  take  one  of  the  two  MSK  frequencies  as 
carrier,  and  study  the  resulting  modulation  function.  It  proves  useful  to 
take  the  upper  frequency,  uiQ  +  v,  as  the  alternative  carrier,  for  which  the 
modulation  function  will  be  called  Z'(t).  The  corresponding  signal  modulation 
is 

S'(t)  =  e“iVtS(t) 


since,  obviously 


S'(t)e1(u)o  +  V)t  =  S(t)eiV 


Substituting,  we  find 


S'(t)  =  e"iVt  E  ina  P(t-nA) 
n 


_  -iv(t-nA)„/t  .. 

=  E  a  e  P(t-nA)  , 

n 


or 


where 


S'(t)  =  E  anP'(t-nA)  , 
n 

P’(t)  =  e_iVtP(t) 

-A  <_  t  <_  A 
otherwise 


|d  H-  e-2iVt) 
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It  is  easy  to  see  how  MSK  results  from  this  modulation.  At  any  time, 
only  two  sucessive  terms  contribute  to  S'(t),  and  if  the  corresponding  code 
bits  are  equal,  the  upper  frequency  results,  since 

P'(t)  +  P'(t-A)  =  1. 

If  the  bits  are  different,  the  lower  frequency  results  because 
P'(t)  -  P'(t-A)  =  e_2iVt, 

and  this  modulation  converts  the  carrier,  tOQ  +  v,  into  the  lower  frequency, 

CO  -  V 

o 

The  representation  using  S'(t)  also  explains  a  common  method  of  generation 
of  MSK.  It  can  be  directly  verified  that 


P'(t)  =  Jn(s)  PQ(t-s)  ds  , 


where 


H(s)  =  vPQ(s)e 


The  proof  follows  easily  from  the  expression 
A/2 


P'  (t)  =  V  |  e 
-A/2 


1 


<-♦4 


-2ivs 


P  (t-s)ds  =  ve 
o 


/ 


-2ivt  /  2i  vo 


*-4 


P  (0)d0  . 

o 


Thus  P'(t)  is  the  response  of  a  filter  to  the  basic  pulse  PQ(t);  the  impulse 
response  of  the  filter  is  H(s)  (excluding  the  realizability  delay),  which 
represents  a  filter  matched  to  a  burst,  one  chj.p  in  duration,  of  the  lower 
frequency. 
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Using  this  fact,  we  can  write 


S' (t)  -  j H(s) 

which  shows  that  S'(t)  can  be  formed  by  generating  a  BPSK  waveform,  using  the 
upper  frequency  and  code  sequence  a^,  and  passing  this  waveform  through  the 
filter  described  by  the  modulation  function  H(s). 

Both  representations  of  MSK  are  useful,  and  they  suggest  different  reali¬ 
zations  of  the  matched  filter  processor.  In  either  case,  the  filter  output 
SNR  is  determined  by  the  magnitude  of  the  signal  autocorrelation  function. 
Using  the  S-representation  (i.e.  the  carrier  u)q),  we  have 


£  a  P  (t-nA-s) 
n  o 
n 


ds 


C(t) 


/ 


S*(o)S( 


0  +  t)d  o 


=  £  im  n  a  a  [P(c  -  nA)P(o  +  t  -  mA)do 
n  m  | 


=  A  £  iZ  C - p ( t  -  i  A)  . 

In  this  formula,  is  the  code  autocorrelation,  as  before,  and  p(t)  is  the 
normalized  autocorrelation  function  of  the  MSK  pulse: 

p(t)  =  iJp(0  )p(0  +  +)d0 


(1  -  2^ ) cos  vt  +  ^  sin(v|t|)  ;  -2A  <  t  <  2  A 


0 


otherwise. 


Note  that  this  function  is  non-zero  for  a  time  four  chips  in  duration.  Thus, 
in  general,  many  terms  contribute  to  the  value  of  C(t)  at  any  given  time,  but 
for  t  =  0, 


C(0)  =  A  CQp(0)  =  L  A  , 


since  p(£A)  vanishes  for  1,  and  the  terms  £  =  1  and  £  =  -1  cancel  (C^  and 

p(£A)  are  even  functions  of  £)  .  Had  we  used  the  S'-representation  (carrier 

a)  +  v) .  we  would  have  found  the  autocorrelation  function 
o 

C'(t)  =  e"1VtC(t)  , 
but  in  either  case  the  output  SNR  is 


SNR(t) 


ES  I  C ( t  -  T)  I 

N  I  C(0) 
o 


for  a  signal  with  arrival  time  T  . 

The  matched  filter  structure  follows  exactly  as  in  the  BPSK  case.  Start¬ 
ing  with  the  S-representation,  we  define  the  prefiltered  modulation  function 
ZLCt) : 


zx(t) 


Jn 


a  )Z(  a  +  t)do 


and  find  that 


W(t)  =  f 


* 

S 


(  a  )Z(  a  +  t)da 


=  Eai  Z,(t  +  n  A) 
n  1 
n 
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Note  that  Z^(t)  is  formed  by  passing  Z^(t)  through  a  filter  matched  to  P(t), 
considering  coo  to  be  the  carrier.  This  can  be  realized  as  a  bandpass  filter 
or,  alternatively,  the  filtering  could  be  performed  on  baseband  signals  ob¬ 
tained  using  the  average  frequency  as  the  reference.  A  problem  arises  in  the 
next  step,  the  tapped  delay  line  structure,  since  the  weights,  a^i  n,  are 
complex.  To  implement  this  literally  at  baseband  requires  a  pair  of  tapped 
delay  lines,  in  which  the  even-numbered  taps  from  one  line  are  weighted  and 
combined  with  the  weighted,  odd-numbered  taps  of  the  other  line  to  form  one 
output.  The  remaining  taps  of  the  two  lines  are  used  to  form  the  other  output, 
representing  the  real  and  imaginary  parts  of  W(t),  according  to  the  equations 


Re{W(t) }  =  Z  (-l)na2n  X1(t  +  2n  A)  +  K-l^a^^  [t  +  (2n+l)A] 

n  n 

ImfW(t)}  =  T.  (-Dna2n  Yx(t  +  2^  A)  -  Z(-l)na2  Xj.  [t  +  (2n+l)A], 
n  n 

where  and  Y^  are  the  I  and  Q  components  of  Z^ : 

Zx(t)  =  Xx(t)  +  i  Yx(t)  . 

The  same  effect  has  been  very  closely  approximated  in  a  SAW  structure, - 
operating  as  a  bandpass  filter  and  delay  structure,  in  which  each  element  of 
delay  is  decreased  by  one-fourth  the  local  carrier  wavelength  to  effect  the 
90-degree  phase  shifts  required  by  the  successive  factors  i  n  . 

Another  matched  filter  structure  results  from  using  the  S-representation, 
where  the  received  modulation  is  Z' (t)  and  the  upper  MSK  tone  is  considered 
to  be  the  carrier.  We  immediately  find  that 


W'  (t)  =  J S'*(o  )  Z' (O  +  t)do 

=  Z  anZ^(t  +  n  A)  , 
n 
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where 


Z|(t)  =  y"p'*(a  )  Z'(a  +  t)do  . 

The  delay  line  structure  is  now  identical  to  that  for  BPSK,  with  real  weights, 

a  ,  and  the  complexity  has  been  shifted  to  the  prefiltering  which  produces 
a  ^ 

Z'  (t) .  Since  P'  ( a )  =  P '  (- a  )  ,  we  have 

z[(t)  =  fp'Co  )Z'  (t-o  )do  , 


hence  Zj  (t)  is  obtained  by  passing  Z'(t)  through  a  filter  matched  to  P'(t). 

The  impulse  response  of  this  filter  is  a  sum  of  the  two  MSK  frequencies. 

It  is  easier,  however,  to  represent  this  filter  as  a  cascade  of  two  fil¬ 
ters,  using  our  previous  result 

P'(t)  =  /H(a)Po(t  -  s)ds 

Thus  Z'(t)  can  be  passed  first  through  the  filter  described  by  H(t),  and  then 
through  the  filter  matched  to  PQ(t),  (or  vise  versa)  according  to 

Z|(t)  =  J Jpq(s)  H(a)Z'(t  -  s  -  a)  dods  . 

All  this  filtering  can  be  accomplished  on  a  carrier,  using  appropriate  bandpass 
filters. 

Since  the  real  signal  described  by  modulation  Z|(t)  and  carrier  coq  +  v  is 
identical  to  the  signal  described  by  Z^(t)  and  carrier  0)Q,  it  is  obvious  that 
the  filtering  operations  described  by  P(t)  in  the  S-representation  and  P'(t) 
in  the  S-representation  are  identical  if  carried  out  at  bandpass.  Having  per¬ 
formed  this  bandpass  filtering,  we  are  then  led  to  the  requirement  of  complex 
weights  if  we  use  a)  as  I  and  Q  reference.  The  second  approach  shows  that  this 
problem  simply  goes  away  if  the  upper  tone  is  used  as  reference  in  going  to 
baseband . 
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The  two  representations  also  lead  to  different  structures  If  the  filter¬ 
ing  is  all  done  at  baseband.  Since  P(t)  is  real,  the  components  of  Z^(t)  are 
each  obtained  from  the  corresponding  component  of  Z(t),  using  identical  filters 
as  noted  above.  However,  if  the  incident  signal  is  converted  directly  to  base¬ 
band  using  the  frequency  <d  +  v,  the  resulting  signals  X'  and  Y',  are  the  com¬ 
ponents  of  Z '  : 

Z'(t)  =  X'(t)  +  i  Y ' ( t)  . 

From  these  we  must  obtain  the  components  of  Z|(t),  according  to  the  formulas 
already  given.  If  we  define 

Z'a  (t)  =  Jh(0)  Z'(t-  o)do  , 

then 

Z1  (t)  =  f Po(°  )Za  (t  _0  )d°  • 

Since  PQ(t)  is  real,  each  component  of  Zj  (t)  is  obtained  by  passing  the 
corresponding  component  of  Z'  (t)  through  a  filter  matched  to  P  (t)  . 

3  O 

To  obtain  Z'(t)  requires  four  filters.  In  terms  of  the  I  and  Q  components: 

3 


z'  (t) 

=  X'  (t) 

+  i  Y'(t) 

a 

a 

a 

Z'  (t) 

5  X'  (t) 

+  i  Y'  (t)  , 

we  have 


x;(t) 

Ya(t) 


Jh^C  a  )X'  (t-  a  )dc  -  jH2(o)Y'(t-a)da  , 
f H  (a  )X'(t-a  )d a  +  )Y'(t  -o  )do  . 
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The  filter  impulse  responses  are  obtained  from  H(o) 


H^o  )  = 

h2(o  )  = 


v  Pq ( cr  )  cos (2vo) 

-v  P  ( o  )  sin(2vo  ) 
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V.  SAMPLED-DATA  APPROXIMATIONS  TO  THE  MATCHED  FILTER 


Consider  the  expression 

W(t)  =  E  an  Z  (t  +  nA) 
n 

for  the  matched  filter  output,  in  terms  of  the  prefiltered  modulation,  Z^(t). 
This  was  derived  first  for  BPSK,  and  also  for  MSK,  when  an  appropriate  choice 
was  made  for  carrier  frequency.  Suppose  this  output  is  sampled  once  every  6 
seconds.  The  output  samples  are  then 

Wk  =  W(k6)  =  E  anZ1(k6  +  nA)  . 

n 

Obviously,  if  6  =  A/M,  for  some  integer  M,  then  the  output  samples  depend  only 
on  values  of  Z^(t)  at  multiples  of  6,  and  the  same  output  sequence  can  there¬ 
fore  be  obtained  by  sampling  the  components  of  Z^t)  at  the  rate  1/6,  and  then 
implementing  the  sums  to  form  the  components  of  the  by  discrete  delay  struc¬ 
tures.  The  continuous  delay  line  is  replaced  by  a  sampled-data  tapped  line, 
such  as  a  CCD  device  or  digital  shift  register.  When  provided  with  an  output 
formed  as  a  weighted  sum  of  tapped  signals,  the  structure  is  a  CCD  or  digital 
correlator.  In  this  report  we  discuss  only  the  approximation  of  matched  filters 
by  such  structures,  limited  to  the  use  of  binary  weights.  Note  that  such  a 
structure  is  a  cascade  of  delay  cells,  one  for  each  sample,  and  NM  in  number. 
Only  every  M*”*1  cell  is  tapped  to  form  the  output. 

Since  only  sampled  outputs  are  available,  there  is  no  guarantee  that  a 
signal  component  will  be  sampled  at  the  moment  of  peak  output  SNR.  This  is 
the  well-known  sampling  loss,  which  we  characterize  by  writing  the  SNR  of 
sample  as 

SNR(W,  )  =  (E  /N  ) .  L.  . 
k  s  o  k 

The  largest  value  of  L^,  which  represents  the  implementation  loss  of  the 
sampled-data  approximation,  will  be  a  function  of  "sampling  phase”,  the  inter¬ 
val  between  the  actual  signal  arrival  time  and  the  nearest  sampling  time. 
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For  either  BPSK  or  MSK,  we  know  that  the  output  SNR  of  the  matched  filter  is 
SNR(t)  -  «,/»„).  |-%^-|  . 

where  t  is  the  signal  arrival  time,  and  C(t)  is  the  signal  autocorrelation 
function.  If  the  sampled-data  receiver  uses  the  same  prefilter  to  form  Z.(t), 
then 


L 


k 


C(k6  -  t | 2 
C(0)  1 


Suppose  the  nearest  sampling  time  to  the  time  of  signal  arrival  is  kQ6,  and 
that 

k  6  =  T  +  s 
o 


with 


<  s  < 


5 

2 


Then 


k  +j 
o  J 


C(j 6  +  s) 


C(0) 


and  the  peak  value  will  occur  for  j  =  0: 


Lk  (s) 
o 


C  (s) 
C(0) 


sj  <  6/2 


If  we  neglect  the  code  autocorrelation  sidelobes  which  enter  into  C(s) 
for  | s | <  6/2  (these  are  only  the  near-in  sidelobes,  which  will  be  small  for 
codes  of  practical  interest),  we  obtain 


Lk  (S^  =  P02^ 

o 

for  BPSK,  and 

(s)  =  p2(s) 

o 
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for  MSK.  Obviously,  sampling  loss  can  be  reduced  by  sampling  at  a  high  mul¬ 
tiple  of  the  chipping  rate,  but  it  is  desirable  to  find  other  ways  of  control¬ 
ling  this  loss,  so  that  sampling  rates  can  be  as  low  as  possible  (or  chipping 
rates  as  high  as  possible  for  a  given  sampling  rate). 

Another  way  to  control  sampling  loss  is  to  depart  from  the  prefilter  asso¬ 
ciated  with  the  original  matched  filter  receiver.  This  is  reasonable,  since 
the  discrete  delay  structure  approximates  the  continuous  version  of  the  matched 
filter,  and  it  is  no  longer  clear  that  the  original  prefilter  is  still  optimum. 
Nothing  can  be  lost,  in  any  case,  by  reopening  the  issue  of  optimization,  but 
we  must  now  recognize  that  sampling  phase  enters  as  a  parameter  of  performance. 

Before  discussing  optimization  criteria,  we  characterize  the  class  of  re¬ 
ceiver  structures  to  be  considered.  These  represent  a  simple  generalization 
in  which  the  "prefiltered  modulation",  Z^(t),  is  defined  by 

Z^t)  =  /G*(c)Z(0+t)dc, 

which  represents  the  output  of  some  filter,  described  by  G(o  ),  to  the  input 
signal.  The  output  sequence  is  defined  as  before,  by 

W  =  Z  a  Z,  (k.5  +  n  A)  , 
k  n  1 

n 

and  hence  the  structure  is  appropriate  to  BPSK,  with  G(a  )  replacing  PQ(o  ), 
and  also  to  MSK  in  the  S ' -representation ,  with  G( o )  replacing  P'(o  ).  Since 
performance  is  characterized  by  SNR  alone,  in  this  analysis,  we  first  compute 
output  SNR  for  the  generalized  receiver,  beginning  with  BPSK. 

Expressing  the  input  modulation  as  a  sum  of  signal  plus  noise,  we  have 

Z^t)  =  Ae1®  L  an  g(t  -  T  -  nA)  +  N^t  , 

n 

where 

g(t)  =  /g*(o)Po(o+  t)do 
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and 


N^t)  =  /c*(  0  )N(  o  +  t)d  i 


It  follows  that  the  signal  component  of  W^  is 

Ae16  E  Cff  g(k6  -  HA  -  t) 
i  * 

We  neglect  the  contribution  of  code  sidelobes  to  this  term,  and  put  t  =  k  6-s, 

as  before,  so  that  the  signal  component  of  W,  is 

o  3 

LAe1^  g( j 6  +  s)  . 

The  covariance  of  the  noise  process,  N^(t),  is  easily  computed: 

E  N*  (t)N1(t')  =  2N  ,//«  o)G*(o')6(o+  t  -  a'-t')dcdo' 

=  2Nq  Jg*  (o'  )G(a'  +  t*  -t)dcf' 

and  hence  the  noise  component  of  any  output  sample  has  mean  squared  value 

2NqL  J  |g( o ) |2  d  a 

when  code  sidelobes  are  neglected  as  before.  Combining  these  expressions 
yields  the  desired  signal-to-noise  ratio  loss 

+  s)  |2 


A/iG(or 


Our  objective  now  is  to  optimize  performance  by  appropriately  choosing 
the  function  G( a  )  which  describes  the  prefiltering.  We  are  free  to  direct 
the  optimization  to  the  output  sample  kQ  (i.e.  put  j=0) ,  since  realizability 
can  always  be  assured  by  the  later  addition  of  a  delay.  Thus  G( o )  will  be 
chosen  to  optimize 
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] 


Lk (s) 


| Jg* (o  )Pq(  a  +  s)do|2 
&J\G(o  )  1 2d  a 


If  the  sampling  phase,  s,  were  fixed  and  known,  then  the  optimum  choice  would 
be  G(o  )  =  P  ( o +  s),  by  an  obvious  application  of  the  Schwarz  Inequality. 

This  is  again  the  matched  filter  of  course,  with  just  the  right  delay  included 
in  the  prefilter  to  assure  that  the  k  ^  sample  occurs  at  a  time  of  peak  output. 

The  simplest  approach  to  the  problem  of  optimizing  L,  (s)  is  to  think  of 

ko 

the  sampling  phase,  s,  as  a  random  variable,  and  to  replace  the  signal  component 
of  Z (t)  by  its  average  over  sampling  phase.  This  has  the  effect  of  replacing 
the  factor  P  (o  +  s)  by  the  average, 


6/2 


PQ(0)  5  (1/6) 


nf?0(a  + 

-6/2 


s)ds 


in  the  formula  for  SNR.  The  corresponding  loss  factor  is  then 


|G*(0  )Pq(o  )d  o|  2 
Aj"|G(o  )  j2do 


This  factor  is  maximized,  again  by  the  Schwarz  Inequality,  by  the  choice 


G  (  a  )  =  Pq  ( a  ) 

which  implies  a  prefilter  matched  to  the  "average  pulse",  P0(°  )•  The  actual 
loss  function  of  a  receiver  using  this  prefilter  is,  of  course, 


Lk  (S> 


|/pq(  o  )Po(  o+  s)da|2 
a/[Po(°  )]2d0 
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The  impulse  response  of  this  filter  is  easily  obtained: 


o  +  | 


C(o)  =  (1/6) 


/ 


6  p0(s)ds 
0  "  2 


o  < 


a  -  a 


-  2 


<  I  o  I  <  A  +  6 


-  2 


A  +  6  | 

— ^ —  <  o 


If  6  is  very  small,  i.e.  ,  many  samples  per  chip,  the  portion  where  G  (o)  varies 
linearly  becomes  very  small  and  G( a  }  approaches  P  (o),  leading  back  to  the 
matched  filter.  At  the  other  extreme,  when  6  =  A,  G(o  )  is  a  symmetrical 
triangular  function,  2A  wide  at  the  base. 

The  loss  function  for  this  filter  is  easily  computed,  either  by  direct 
evaluation  or  by  use  of  the  relation 


/ P  (o  )P  (cr+s)do  = 
Jo  o 


s+  | 


if  po(0)do 


Strictly  speaking,  (s)  was  defined  only  for  values  of  s  in  the  range 
| s | _<  6/2,  corresponding  to  SNR  of  the  output  sample  with  maximum  signal.  It 
should  be  obvious  from  our  definition,  however,  that  s  can  take  any  value  and 
the  loss  (s)  will  describe  the  SNR  of  a  particular  output  sample.  For  ex¬ 
ample,  if  6/2  <  s  <  3  6/2  , 


\  <s>  =  Lk0+l(s  "  °  ) 


which  characterizes  the  SNR  of  output  sample  W,  +  j  etc 
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The  curves  of  Fig.  V-l  show  the  variation  of  SNR  with  sampling  phase 
for  this  filter,  for  the  two  cases  M  =  1  (6  =  A)  and  M  =  2  (6  =  A/2).  The 
curve  corresponding  to  the  original  matched  filter  is  also  shown.  Expanded 
plots  of  SNR  loss,  in  dB,  are  presented  in  Fig.  V-2.  Explicit  formulas  are 
obtained  in  Appendix  A-I. 

A  second  approach  to  optimization  is  to  choose  G( a )  to  maximize  the  aver¬ 
age  SNR,  averaged  over  sampling  phase,  i.e. 

5/2 

°  -6/2  ° 

This  is  a  somewhat  more  difficult  problem,  but  an  explicit  solution  can  be 
found,  and  a  derivation  is  presented  in  Appendix  2.  The  corresponding  filter 
is  almost  identical  to  the  one  just  derived,  with  almost  identical  performance 
Average  SNR  losses  (in  dB)  are  given  in  Table  I  for  three  filters: 


"A"  Filter  maximizing  average  SNR  » 

"B"  Filter  matched  to  average  pulse  PQ(t)  ,  and 
"C"  Filter  matched  to  the  pulse  Pq  (t). 


TABLE  I 


AVERAGE  SNR  LOSS  IN  dB  FOR  THREE  FILTERS 


M 


"C 


II 


1 

2 

4 


1.70 

0.78 

0.38 


1.82 

0.80 

0.38 


2.67 

1.19 

0.57 


Worst-case  losses  for  filters  "B"  and  "C"  are  as  follows: 

TABLE  II 

WORST-CASE  SNR  LOSSES  (IN  dB) 

M  "B"  "C" 

1 
2 
4 


4.26 

1.71 

0.78 


6.02 

2.50 

1.16 


TR-550  (V-1) 


Fig.  V-1.  BPSK  Loss  Factor  for  three  filters 


Fig.  V-2.  BPSK  Sampling  Loss  in  dB. 
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It  is  also  possible  to  design  a  filter  which  yields  a  loss  independent 
of  sampling  phase,  for  the  one  output  sample  having  the  largest  SNR.  The 
penalty  paid  for  this  uniformity  is  a  relatively  large  loss,  as  follows: 

TABLE  III 

SNR  LOSSES  FOR  FILTER  HAVING  CONSTANT  LOSS  (IN  dB) 

M  Constant  Loss 

1  3.01 

2  1.25 

4  0.58 

Derivation  of  this  filter  and  its  performance  is  presented  in  Appendix  3. 

A  broader  class  of  sampled-data  receivers  emerges  if  every  cell  of  the 
delay  structure  is  tapped,  forming  an  output  of  the  form 


Wk  =  E  A£  Z1[(k+R)6] 

The  sum  would  range  over  NM  terms,  corresponding  to  the  signal  duration.  A 
sub-class  results  if  the  weights,  A^  ,  are  equal  in  groups  of  M,  and  corre¬ 
spond  to  the  code  bits,  so  that 

M 

WR  =  1  an  1  ZjdcS  +  nA  +  m<S) 
n  m=  1 


This  output  can  be  written 


where 


W.  =  Y  a  Z-(k<5  +  nA) 
k  n  2 

n 


M 

Z2( kA)  =  T.  Zj  ( (k  +  m)6] 
m«l 


In  other  words,  the  summation  over  the  M  terms  per  chip  can  be  carried  out 
once  and  for  all,  ahead  of  the  delav  structure,  since  these  terms  always  have 
equal  weight.  Moreover,  the  same  result  can  be  achieved  ahead  of  the  sampling 
step  by  forming 
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M 

Z0(t)  =  ?:  Z  (t  +  m  &) 
m=  1  1 

and  sampling  Z.,(t).  Finally,  substituting  for  Zj(t),  we  obtain 

M  f  * 

Z.,(t)  =  T.  I  C,  (  o  )  Z  ( O  +  t  +  m  '  )do 
m=  L  * 

=  Jc*  (  o)  Z  (  0  +  t )  d  0 
M 

where  l!  ,  (  0  )  V  C.(o  -  m  <:)  . 

m=  1 

In  other  words,  the  effect  of  the  multiple,  equal-weight  taps  can  be 
achieved  in  the  prefiltering,  and  hence  this  subclass  offers  no  performance 
advantage,  since  we  have  already  optimized  the  prefilter.  It  is  interesting 
to  note  that  the  til'  r  we  obtained  by  matching  to  the  average  pulse,  F^ (  c  ), 
can  be  written  in  the  form  of  (;,(o  )  above,  namelv  a  sum  of  M  terms,  each  a 
displaced  version  of  the  simple  triangle  function: 


0 


othe rw i se 


>  i  iusive  of  a  i  ixed  delay).  Thus  a  receiver  using  a  prefilter  characterized 
h ■  tuis  simple  triangular  impulse  response,  in  conjuction  with  a  delay  struc- 
t  .to  tap; "ing  everv  cell  in  the  special  wav  described  above,  would  be  exactly 
c  ;  .i  valent  to  tin*  receiver  we  discussed  earlier,  with  prefilter  "B". 

!  he  same  approach  to  sampled-data  receiver  optimization  can  be  applied 
t  -  Msg,  beginning  with  the  generalized  prefiltered  modulation 

/-[  (t)  =  Jc*(  C  )  7.'  (  c  +  r)  d  ■  . 

:,av«-  ret  lined  the  prime  in  our  notation  to  emphasize  the  assumption  of  the 
uigh.  r  rrequenev,  .  +  v,  as  carrier  frequency.  The  remaining  part  of  the 

;  ■  t  '.easing  ;  s 
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W!  =  E  a  7' (k6  +  nA) 
k  n  1 

n 

obtained  as  before  by  sampling  the  continuous  output  at  times  k.6.  The  calcu¬ 
lation  of  output  SNR  is  exactly  parallel  to  the  BPSK  case  just  treated,  with 
the  resulting  loss  function 


|  Jg  (  a  )P*  (  o  +  s)d  o  |2 

aJ|G(  a  )  |2  da 


As  before,  the  simplest  optimization  procedure  is  to  match  to  the  "average 
pulse",  but  we  have  a  choice  now  of  doing  this  in  either  of  the  two  representa¬ 
tions,  and  the  results  are  not  quite  the  same.  Averaging  in  the  S' -representa¬ 
tion  means  choosing 


o  +  6/2 


G(  a  )  =  P 


'  (  o  )  =  J  P'  (u)du 


c  -  6/2 


while  averaging  in  the  S-representation  corresponds  to  the  choice 


G(  0  )  =  e’iVO  P  (  o  ) 


c  +  6/2 


iV°  ’  ^ j 


c  -  6/2 


Al  though 


_ .  ,  ,  - ivu  „ ,  . 

P'  (u)  =  e  P(u)  , 


the  resulting  filters  are  different,  and  it  happens  that  the  second  choice 
(S-representation)  yields  better  performance.  The  performance  analysis  is 
given  in  Appendix  A,  and  the  loss  curve  for  M  =  1  is  shown  in  Fig.  V-3  compared 
to  the  result  for  the  choice:  G(o  )  =  P'  ( a  ) ,  which  is  matched  to  the  MSK  pulse 
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itself.  It  is  seen  that  the  results  are  very  close,  and  similar  curves  for 
M  >  2  show  insignificant  differences  between  matching  to  the  MSK  pulse  and 
matching  to  an  averaged  MSK  pulse.  Average  losses  are  shown  in  Table  IV  for 
filter  "B":  matched  to  the  average  pulse,  and  filter  "C":  matched  to  the  MSK 
pulse  itself: 


TABLE  IV 

AVERAGE  SNR  LOSS  FOR  TWO  MSK  FILTERS  (IN  dB) 

M  "B"  "C" 

1  0.742  0.825 

2  0.204  0.213 

4  0.054  0.055 

No  explicit  solutions  have  been  obtained  to  the  problems  of  optimizing 
average  SNR  and  of  obtaining  a  constant  loss  function,  in  the  MSK  case.  In 
view  of  the  small  actual  improvement  of  filter  "B"  above,  relative  to  the 
conventional  choice  "C",  there  is  probably  little  to  be  gained  by  these  other 
approaches  to  optimization. 

By  the  very  nature  of  its  autocorrelation  function,  p(t),  MSK  shows  less 
sensitivity  to  sampling  phase  than  does  BPSK.  There  is  little  to  gain  by 
filter  optimization,  and  only  about  0.5  dB  average  improvement  of  double-rate 
(M  =  2)  sampling  over  the  single-rate  (M  =  1)  case.  This  compares  with  a  1  dB 
to  1.5  dB  improvement  (depending  on  filter  choice)  in  the  BPSK  average  SNR, 
when  going  from  single-rate  to  double-rate  sampling.  Worst-case  SNR  losses 
are  also  milder  for  MSK.  The  following  Table  compares  the  worst  case  losses 
of  the  filter  "B"  for  BPSK  (the  "recommended  filter",  matched  to  an  average 
BPSK  pulse)  with  "C"  for  MSK  (the  conventional  MSK  filter,  matched  to  the 
MSK  pulse  itself)  : 

TABLE  V 

WORST-CASE  SNR  LOSSES  (IN  dB) 

M  BPSK  "B"  MSK  "C" 

1  4.26  2.44 

2  1.71  0.63 

4  0.78  0.16 
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We  have  found  that  the  sampled-data  MSK  receiver  uses  the  same  discrete 
delay  structure  as  BPSK,  where  each  cell  delays  by  one  sample,  and  the  struc¬ 
ture  has  one  tap  per  chip.  For  BPSK  we  showed  that  every  cell  could  be  tapped 
with  weights  equal  in  groups  of  M,  and  the  result  was  the  same,  provided  a 
suitable  change  was  made  in  the  prefilter.  Specifically,  the  receiver  with 
prefilter  g(o),  and  M  taps  per  chip  was  equivalent  to  a  receiver  having  one 
tap  per  chip  and  prefilter 

M 

G  ( a  )  =  Z  g  (  a  -  m<5)  ,  <5  =  A/M  . 

m=l 

The  same  result  holds  for  MSK,  since  the  original  proof  depended  only  on  the 
architecture  of  the  receiver,  and  not  on  the  signal  modulation  itself.  For  a 
given  receiver  tapping  every  sample,  it  is  easy  to  find  the  equivalent  pre¬ 
filter  corresponding  to  one  tap  per  chip,  since  we  start  with  g(  O )  and  find 
G(  c )  by  the  above  formula.  It  is  not  always  easy  to  go  the  other  way  however 
representing  a  given  G(  a  )  in  the  form  required  to  find  g  (a).  It  is  per  laps 
of  interest  to  note  that  this  problem  can  be  solved  for  the  "C"  filter,  since 

G( 0  )  =  P'  ( a  )  =  j H(s)  Pq(s  -  u)  du. 

We  simply  represent  the  basic  pulse,  Pq,  as  a  sum  of  M  pulses,  back-to-back, 
each  lasting  one  sample  instead  of  one  chip,  and  substitute  in  this  integral. 
The  result  is  a  representation  of  G(o  )  in  the  desired  form.  The  receiver 
prefilter  was  factured  into  a  cascade  of  two  filters,  with  impulse  responses 
H ( s )  and  P  (s) ,  hence  it  is  only  necessary  to  change  the  second  from  a  filter 
matched  to  a  pulse  one  chip  in  duration  to  a  filter  matched  to  a  pulse  one 
sample  in  duration.  Combined  with  a  change  from  one  tap  per  chip  to  one  tap 
per  cell,  the  result  is  the  equivalent  receiver. 
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VI.  CONCLUSIONS 

The  two  modulations  discussed  in  this  report,  BPSK  and  MSK,  are  common 
examples  of  spreading  waveforms  useful  in  a  wide  range  of  spread-spectrum 
systems.  A  great  variety  of  signals  can  be  designed,  using  these  coded  pulses 
as  basic  building  blocks,  hence  the  results  of  this  study  have  a  wide  range  of 
application.  Also,  it  appears  that  sampled-data  receiver  structures  of  the 
kind  treated  here  are  well  suited  to  these  waveforms,  where  the  major  part  of 
the  processing  gain  is  achieved  by  coherent  integration  of  chips  in  a  CCD  or 
digital  correlator. 

One  result  of  this  study  is  a  systematic  method  for  the  design  of  sampled- 
data  receivers  which  can  easily  be  extended  to  other  spread-spectrum  modulations. 
This  method  begins  with  the  exact,  continuous-time  matched  filter,  and  the  dis¬ 
crete  portion  of  the  receiver  is  then  inferred  from  the  requirement  that  its 
output  samples  should  be  identical  with  samples  of  the  true  matched-filter  out¬ 
put.  The  second  step  of  the  design  process  is  to  choose  the  prefilter,  imple¬ 
mented  at  bandpass  or  base  band  in  analog  form,  to  optimize  SNR  performance  as 
a  function  of  sampling  phase. 

The  other  result  of  this  analysis  is  the  design  and  performance  of  a  num¬ 
ber  of  prefilters  for  use  with  3PSK  and  MSK.  The  differences  in  performance 
are  not  great,  but  the  few  dB  they  offer  can  in  some  cases  be  effectively 
traded  off  against  some  other  parameter  which  may  be  near  a  practical  bound 
in  value.  For  example,  on.  would  ordinarily  not  use  BPSK  with  only  one  sample 
per  chip,  but  the  fact  that  the  6-dB  worst-case  loss  associated  with  the  con¬ 
ventional  filter  can  be  reduced  to  a  3-dB  loss  (independent  of  sampling  phase) 
might  make  this  choice  interesting  in  some  application  in  which  the  sampling 
rate  cannot  easily  be  doubled.  In  another  case,  with  MSK  and  an  unconven¬ 
tional  filter  (filter  "B"),  double-rate  sampling  buys  only  0.5  dB  in  average 
SNR  performance,  compared  to  single-rate  sampling,  a  difference  which  could 
perhaps  be  achieved  more  easily  by  changing  some  link  parameter,  rather  than 
doubling  the  clocking  rate  of  a  substantial  portion  of  the  receiver  hardware. 
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APPENDIX  1 

If  we  define 

PG(S) 


f  Po(s)1  2 

i/[Po(o)i2d0 


t  hen 

Lk  (s)  = 

o 


One  immediately  finds  that 


i/i 


po(0)l' 


'do  = 


1  - 


1_ 

3M 


where  M  =  A/6 


Since  p  (-s)  =  p  (s)  ,  we  evaluate  p  (s)  for  positive  s 
o  o  ° 

vt  have 


When 


<  A  - 


we  have 


When  0 
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APPENDIX  2 

For  an  arbitrary  filter,  G(o),  the  SNR  as  a  function  of  sampling  phase 
is  proportional  to  the  ratio 

\Jg(  o  )  PQ(  O  +  s)d  a  |  2 
f \G(  a  ) | 2d  o 


We  now  wish  to  choose  G(o  )  in  order  to  maximize  the  average  value  of  this 
expression,  treating  sampling  phase  as  a  uniformly  distributed  random  variable. 
The  solution  to  this  problem  will  be  sketched  briefly  here.  The  numerator  is 


// 


C  (o  )C,(  o,)Pq(  o+  s)P  (o'  +  s)dod  o' 


and  its  average  can  be  expressed  as 

ffG*(  0  )G(  d  )A(  a  ,  d  )d  a  d  o' 


where 


A(  0  ,  d  ) 


6/2 


P  (c+s)P  (o'  +  s)  ds 
o  o 


-6/2 


The  problem  of  maximizing  the  ratio 


// 


G  (o)G(o')A(o,o,)dodo' 


f |g( a  ) 


"d  o 


is  a  standard  one,  whose  solution  is  to  take  G( o  )  proportional  to  that  eigen¬ 
function  of  the  operator  A,  corresponding  to  the  largest  eigenvalue  of  A.  The 
eigenfunctions  of  A  are  defined  by  the  equation 
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Af(  o  ) 


Ja(  a  ,  d)t(d)dd  = 

Together,  they  span  a  subspace  of  functions,  and  it  is  clearly  disadvantageous 
to  choose  a  G(  o )  which  has  a  component  outside  this  subspace,  since  this  com¬ 
ponent  increases  the  denominator  of  our  ratio  without  helping  the  numerator. 

This  heuristic  argument  is  easily  made  rigorous,  and  once  confined  to  the  sub¬ 
space  in  question,  G( 0 )  can  be  expanded  in  eigenfunctions: 

G(o)  =  E  8  ip  ( a  ) 
m  m 
m 

The  positive  definite  operator  A  will  turn  out  to  have  a  discrete  spectrum, 
and  we  will  have 

//•*  (  o  )A(  o  ,  d  )G(  o*  )d  o  d  o  ' 

while 

/|G(o)|2do  =  I|ej2 
m 

(we  have  assumed  the  ^  to  be  o.  thogonal  and  normalized),  if  ^  fs  the  iargest 
eigenvalue  (they  are  positive),  then  the  choice  G(o  )  =  fQ(a)  solves  our  problem. 

Since  p0(t)  is  a  simple  pulse,  it  is  easy  to  evaluate  A(  o  ,  d )  explictly, 
although  it  is  complicated  to  write  down.  As  a  function  of  o  ,  A(  o  ,  o')  is 
constant  when  |o  |  <  ^  (A  -  5)  and  zero  when  [  o[ >  y  (A  +  6).  Between  these 
bounds,  A ( a  ,  o' )  is  either  zero,  or  linear  (increasing  or  decreasing)  in  a. 

It  follows  that  the  eigenfunctions  are  zero  when  |  o  | >  4  (A  +  6)  and  constant 

ill  i 

when  |  0  |<  j  (A  -  6).  Only  the  portion  of  the  eigenfunctions  in  the  remaining 

range  needs  to  be  determined,  and  by  differentiating  the  integral  equation 

twice,  a  differential  equation  resuits.  This  equation  is  simply  the  equation 

of  a  sinusoid,  and  by  matching  boundary  conditions,  a  discrete  spectrum  is 

obtained.  The  solution  with  maximum  eigenvalue  joins  the  constant  portion  to 

the  value  zero  at  the  ends  of  the  range  with  a  portion  of  a  sinusoid  that  differs 

only  very  slightly  from  a  straight  line,  the  corresponding  solution,  matched 

to  the  averaged  pulse.  The  SNR  attained  is  determined  by  the  eigenvalue  itself, 

whose  evaluation  leads  to  the  values  given  in  the  text. 
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APPENDIX  3 

To  obtain  a  SNR  independent  of  sampling  phase,  it  is  necessary  to  choose 
G(  O )  so  that 


Jg (o)  PQ(o+s)d  a 
G(o)do 

+  s 

is  independent  of  s,  for  |s  |£  6/2  .  By  differentiating  this  expression, 

it  is  seen  that  we  require 


2 


G(j+s)  =  G(-  j  +  s)  ,  -|lsl| 

Thus,  we  can  put 

G(  G  )  =  a(  o  )  for  |  o  |  <_ 


and 


G(f  + 


o)  «  G(-  f  + 


a)  =  b ( a  ) 


for  |  o  |<_  6/2  , 


where  a(  0  )  and  b( o )  are  arbitrary,  and  G( a )  is  taken  equal  to  zero  elsewhere. 
Since  SNR  is  now  independent  of  s,  we  can  choose  s  =  -6/2,  and  the  loss  factor 
will  be 


L 


+ 


2 
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6/2 


A  -  6 


/ 


a  (o)d  o  +  2 

A  -  6  -6/2 


jb2(o)  do 


It  is  not  hard  to  show,  by  means  of  the  calculus  of  variations, 
ratio  is  maximized  when  a(  o  )  and  b(o)  are  constant,  with  b  =  a/2, 
the  filter  referred  to  in  the  text,  and  the  loss  factor  resulting  is 


that  this 
This  is 
given  by 
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APPENDIX  4 

The  basic  expression  for  the  loss  factor  can  be  written  in  the  form 


<u<s>  ■ 

IVs)l2 

°M 

where 

Vs)  5 

I Jg  *(  o  )P(o  +  s)do 

and 

)|2do 

6 

—  l  V 

<*<°>  * 

P(a  )  =  -^  /  P(  o+  u)du 

2 


The  MSK  autocorrelation  function  is 


P(t)  E  +  t)do 


V2A7  77 


'/■i  tk  .7Tt.  ,  1  ,  .Til  1 1  I  I 

(1  -  cos  <^)  +  -  Sin  (-2X  ;  |t|  <  2A 


;  otherwise 


In  terms  of  p(t),  we  have 


y°>  ■  if 


V 


6 

2 

/ 


P (  a  +  u)du 


P(  o  +  s)da 


s  +  | 


P(s-u)du  *  -r 
0 


i/ 


p(u)du 


8  "  2 


42 


and  also 


—  ^  0  +  U)P(°+  u')dudu* 


d  o 


~2 jj P(u'  "  u)dudu'  =  ~2  / 


p(w) (6  -  w)dw 


In  this  last  integral  we  put  w  =  xo  and  substitute  for  p(w) ,  with  the  result 


°M  =  2 


(!  -  fj^cosC^)  +  ?  sin(§>  1  d  ~  x>d*> 


where  M  *  A/6  .  Direct  evaluation  yields 


D  =  ^ 

M  2 

IT 


2 (M  +  1)-(2M-1)cos(2^)-^  sinC^) 


For  M  =  1  ,  this  simplifies  to 


“x  - 

7T 

Turning  to  B^Cs) ,  we  note  that  it  is  an  even  function  of  s  and  that 


F(u)  =  f/ju-  fA)cos(ir)+  i  sin(§r: 


du' 


can  be  evaluated  as 
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* 


1  U,  ,  .TTu  .  4M  .TTu  . 

F(u)  =  7I  (2M-  6)  sin(2M6)  2  cos<2m6  ' 

TT 

When  s  >  2  A  +  ^  ,  we  have  B^s)  =  0,  since  the  averaging  interval  is 

beyond  the  range  of  p(u).  When  2A  -  s  £  2A  +  j  ,  we  have 

Bm(s)  =  F(2A)  -  F(s  -  |  ) 


T  +  sin(^)  j  ^  sin(4 J)-  ^  (2M  +  |  -  -^)cos(^) 

TT  [  TT 


.  TTS ,  I  AM  ,TT  .  1  1  MS,  -  TT  ) 

+  cos(^)  |  cos  (4M>  +  i  (2M  +  2  "  A~  sin  4M  j 


6  & 

Next,  when  -g-  <  s  _<  2A  -  -j  . 


Bm(s)  =  F(s  +  y)  -  F(s  -  ) 


‘  .inCgTj^^n  4)  -i'os  4>} 

(  TT 


,TTS,  )  8M 


,  2M  ,tts,.,„  s,  .  , TT  , 

+  ~ it  COS^2A  ^  “  A  sln(4M  ’ 


Finally,  when  0  <  s  <  |  ,  we  obtain 


Bm(s)  =  F(s  +  |)  +  F(|  -  s>-  2F( 0) 


8M  2Ms  , 

”  2  "  ttA  COS  4M  sln  2A 

TT 


+  14  (^M  -  1) sin  (^)  -  ^  COS  (^)}cos(|f) 


«>)' 


l2A7 
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Hu-  (H-.ik  SNR  always  occurs  in  the  interval 


in  which  case  the  explicit  loss  factor  is 


|”8_M  .  sin(g)cos(4-)+[(4M-l)sin(^j)- 

4M  1 2  ( M  +  1 )  -  (2M-  1  )cos(^j)-  ^sinf^)  | 
For  M  =  1  ,  this  reduces  to 


't 

2  l 


(s)  =  -  — 


+  (r  -1)/2 


COS(27.)"r/2 


sin(E)] 


2  n  -  1 


which  is  shown  in  Fig.  V  -  1. 
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